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Planar brushes of flexible, semiflexible and rodlike diblock polyampholytes are studied using 
molecular dynamics simulations in a wide range of the grafting density. Simulations show linear 
dependence of the average thickness on the grafting density in all cases regardless of different flexi- 
bility of anchored chains and the brushes different equilibrium conformations. Slopes of fltted lines 
to the average thickness of the brushes of semiflexible and rodlike polyampholytes versus the graft- 
ing density are approximately the same and differ considerably from that of the brushes of flexible 
chains. The average thickness of the brush of flexible diblock polyampholytes as a function of the 
grafting density is also obtained using a simple scaling analysis which is in good agreement with our 
simulations. 

PACS numbers: 82.35.Rs, 61.41. -|-e, 87.15.-v 
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I. INTRODUCTION 

Macromolecules containing ionizable groups when dis- 
solve in a polar solvent such as water, dissociate into 
charged macromolecules and counterions (ions of op- 
posite charge). Depending on acidic or basic prop- 
erty of their monomers, ionizable polymers in solu- 
tion can be classified into polyelectrolytes and polyam- 
pholytes. Polyelectrolytes contain a single sign of charged 
monomers and polyampholytes bear charged monomers 
of both signs. These macromolecules are often water- 
soluble and have numerous industrial and medical appli- 
cations. Many biological macromolecules such as DNA, 
RNA, and proteins are charged polymers. In polymer 
science, charged polymers has been an irnportant sub- 
ject during last several decades [E [3, H, Contrary 
to a polyelectrolyte chain in which the intra-chain elec- 
trostatic interactions are repulsive and tend to swell the 
chain, in a polyampholyte chain attractive interactions 
between charged monomers of opposite sign tend to de- 
crease the chain size. Oppositely charged monomers can 
be distributed randomly along a polyampholyte chain 
or charges of one sign can be arranged in long blocks. 
With the same ratio of positively and negatively charged 
monomers (isoelectric condition), behaviors of a single 
polyampholyte and the solution of polyampholytes de- 
pend noticeably on the sequence of charged monomers 
on the chains. For example, it has been shown that 
the sequence of charged amino acids (charge distribu- 
tion) along ionically complementary peptides affect the 
aggregation behavior and self-assembling process in the 
solution of such peptides [1, Also, using Monte Carlo 
simulations it has been shown that charged monomers se- 
quence of charge-symmetric polyampholytes affect their 
adsorption properties to a charged surface 

The properties of the system of polymers anchored on 



a surface are of great interest both in industrial and bio- 
logical applications and academic research. With a suffi- 
ciently strong repulsion between the polymers, the chains 
become stretched and the structure obtained is known as 
a polymer brush. Planar and curved brushes formed by 
grafted homopolymers have extensively been investigated 
by various theoretical methods 8, s|, [13, [lH. The an- 
chored polymers of a brush may be consisting of charged 
monomers. In this case, the brush is known as a poly- 
electrolyte or a polyampholyte brush depending on the 
charged monomers of the chains being of the same sign 
or being composed of both signs. In a brush of charged 
polymers, electrostatic interactions introduce additional 
length scales such as Bjerrum length and Debye screening 
length to the system. 

In a polyelectrolyte brush, the repulsion of electro- 
static origin between the chains can be sufficiently strong 
even at low grafting densities, making it easy for the sys- 
tem to access the brush regime. Polyelectrolyte brushes 
have been investigated extensively using both theoretical 
El [11 [M, 15,^1 eJIlJllJl^ computer simula- 
tion methods Illliallallll ill. At high enough 
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grafting densities and charge fractions of polyelectrolyte 
chains, most of counterions are trapped inside the poly- 
electrolyte brush and competition between osmotic pres- 
sure of the counterions and elasticity of the chains de- 
termines the brush thickness. This regime of a polyelec- 
trolyte brush is known as the osmotic regime in which 
some theoretical scaling methods predict no dependence 
of the brush thickness to the grafting density [H 12^ . 
However, other scaling method that takes into account 
the excluded volume effects and nonlinear elasticity of 
polyelectrolyte chains and is in agreement with exper- 
iment and simulation, predicts a linear dependence of 
the brush thickness on the grafting density [13, [H, [l^ ■ 
Also, it has been shown that diffusion of a fraction of 
counterions outside the polyelectrolyte brush leads to a 
logarithmic dependence of the average brush thickness 
on the grafting density [lB|. 

Electrostatic interactions in a polyelectrolyte brush 
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FIG. 1: (Color online) Sample snapshots of brushes of a) flexible {Ip = 5cr), b) semiflexible {Ip — 25a) and c) rodlike {Ip — 200a) 
diblock polyampholytes at grafting density paa^ — 0.08. Neutral , positively and negatively charged monomers are shown by 
white, red (dark) and blue (black) spheres respectively. Counterions are not shown and periodic boundary condition is removed 
for clarity. 



cause most of the counterions to be trapped inside the 
brush and help the chains to be more stretched and the 
brush to be more ahgned. However, in a brush of over- 
ally neutral polyampholyte chains, most of counterions 
are outside the brush and the electrostatic correlations 
tend to decrease the chains size and the brush thickness. 
At a given value of the grafting density, the average thick- 
ness and equilibrium properties of such a polyampholyte 
brush are mainly determined by the chains properties 
such as fraction and sequence of charged monomers and 
the bending energy. Brushes formed by grafted diblock 
polyampholytes have been investigated by lattice mean 
field modeling [l^, [2^ and computer simulation [23| . The 
effect of chain stiffness, charge density and grafting den- 
sity on spherical brushes of diblock polyampholytes and 
interaction between colloids with grafted diblock polyam- 
pholytes have been studied using Monte Carlo simula- 
tions p7l.l28j. Also, using molecular dynamics (MD) sim- 
ulations, the effects of various parameters such as charged 
monomers sequence, grafting density and salt concentra- 
tion on the average thickness and equilibrium conforma- 
tions of planar semiflexible polyampholyte brushes have 
been investigated poj . 

In this paper, we study planar brushes of flexible, semi- 
flexible and rodlike diblock polyampholytes using MD 
simulations in a wide range of the grafting density. We 
find that in all cases the average brush thickness lin- 
early depends on the grafting density regardless of the 
chains different flexibility. Our results also show that the 
strength of this dependence is considerably weaker in the 
case of the brush of flexible polyampholytes than two 
other cases. Despite mentioned same functionality ob- 
tained for the average thickness versus the grafting den- 
sity for brushes of different polyampholytes we flnd that 



histograms of their equilibrium conformations are notice- 
ably different. Inter-chain correlations are too weak in 
the brush of flexible polyampholytes and the brush prop- 
erties are dominantly determined by single chain behav- 
ior. In this case, dependence of the equilibrium confor- 
mation of the brush on the grafting density is very weak. 
In the cases of the brushes of semiflexible and rodlike 
polyampholytes however, because of the combination of 
electrostatic correlations and strong excluded volume ef- 
fects, collective behavior of the chains is dominant and 
dependence of equilibrium conformations on the grafting 
density is strong. In these cases, we also observe separa- 
tion of the anchored chains into two coexisting fractions. 
Using a simple scaling method which is consistently ap- 
plicable for the brush of flexible chains, we describe theo- 
retically the linear dependence of the brush thickness on 
the grafting density. 

The rest of the paper is organized as follows. In Section 
|TT]we describe our model and simulation method in detail 
and present the results of MD simulations. Our scaling 
analysis to describe dependence of the average thickness 
of flexible diblock polyampholyte brushes on the graft- 
ing density is presented in section IIIII In section IIVI we 
conclude the paper and present a short discussion. 



II. MOLECULAR DYNAMICS SIMULATION OF 
DIBLOCK POLYAMPHOLYTE BRUSHES 

A. The model and the simulation details 

In our simulations which aregerformed with the MD 
simulation package ESPResSo [30| , each brush is modeled 
by M = 25 diblock polyampholyte bead-spring chains of 
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length iV = 24 (24 spherical monomers) which are end- 
grafted onto an uncharged surface at z = 0. The posi- 
tions of anchored monomers which are fixed during the 
simulation, form an square lattice on the grafting surface 

— 1/2 

{x — y plane) with lattice spacing d = pa in which 
Pa is the grafting density of the chains. The fraction 
/ = i of the monomers of each chain are charged and 
the chains consist of an alternating sequence of charged 
and neutral monomers. Each chain contains the same 
number of positively and negatively charged monomers 
with charges e and — e respectively (see Fig. [Ij. Ex- 
cluded volume interaction between particles is modeled 
by a shifted Lennard- Jones potential, 



ULj{r) 



4£ m'' 





i} if r < Tc, 
if r > Tc, 



(1) 



in which e and a are the usual Lennard- Jones parameters 
and the cutoff radius is rc = 2V6(j. 

Successive monomers 
of each chain are bonded to each other by a FENE (finite 
extensible nonlinear elastic) potential [3l|, 



Ubond{r) 







if r > i?o. 



(2) 



with bond strength ktond = 30e/cr^ and maximum bond 
length i?o = 1-5(T. Bending elasticity of the chains is 
modeled by a bond angle potential. 



Ubend{r) = kbend{l ~ COS 6) , 



(3) 



in which 9 is the angle between two successive bond vec- 
tors and kbend is the bending energy of the chains. The 
value of the persistence length, Zp, of the chains depends 
on the value of khend as Ip ~ o. To model brushes 
of flexible, semiflexible and ro dlike chains, we use four 
different values of khcnd namely 0, 5kBT, 25fcBT and 
200kBT respectively. The simulation box is of volume 
L X L X Lz in which L is the box width in x and y direc- 
tions and Lz is its height in z direction and the grafting 
density is given by pa — M/L^. We consider M x N x f 
monovalent counterions to neutralize the chains charge. 
Positive and negative monovalent counterions are mod- 
eled by equal number of spherical Lennard-Jones parti- 
cles of diameter a with charges e and — e respectively. All 
the particles interact repulsively with the grafting surface 
at short distances with the shifted Lennard-Jones poten- 
tial introduced in Eq. [TJ In addition, a similar repulsive 
potential is applied at the top boundary of the simulation 
box and in our simulations Lz = 2Na. All the charged 
particles interact with each other with the Coulomb in- 
teraction 



uc{r) = ksTqiq/— 
r 



(4) 



in which qi and qj are charges of particles i and j in units 
of elementary charge e and r is separation between them. 
The Bjerrum length. Is, which determines the strength of 
the Coulomb interaction relative to the thermal energy. 




FIG. 2: (Color online) The average thickness of the brushes 
formed by flexible {kbend = and SksT), semiflexible 
[kbend = 25kBT) and rodlike {kbend = 200fcsr) diblock 
polyampholytes versus dimensionless grafting density, paO^ ■ 
The lines are linear fits to our simulation data. The slopes 
of solid, dashed, dotted and dash-dotted lines are 41.6, 40.2 
and 22.7 and 20.0 respectively. The size of the symbols cor- 
responds to the size of error bars. 



ksT, is given by Ib — /^ksT, where e is the dielec- 
tric constant of the solvent and we set Ib — 2a 'm our 
simulations. Periodic boundary conditions are applied 
only in two dimensions (a; and y). To calculate Coulomb 
forces and energies, we use the so-called MMM tech- 
nique introduced by Strebel and Sperb [3^ and modified 
for laterally periodic systems {MMM2D) by Arnold and 
Holm 33]. The temperature in our simulations is kept 
fixed at ksT ~ 1.2e using a Langevin thermostat. 

For each value of the bending energy, khend, we do 
simulations of the brush at dimensionless grafting den- 
sities pad^ = 0.02,0.04,0.06,0.08,0.10. In the begin- 
ning of each simulation, all of the chains are straight 
and perpendicular to the grafting surface and all the 
ions are randomly distributed inside the simulation box. 
We equilibrate the system for 1.6 x 10^ MD time steps 
which is enough for all values of the grafting density men- 
tioned above and then calculate thermal averages over 
1500 independent configurations of the system selected 
from 2.25 x 10^ additional MD steps after equilibration. 
MD time step in our simulations is t = O.OItq in which 



To = 



is the MD time scale and m is the mass of 



the particles. 

We calculate the average brush thickness which can 
be measured by taking the first moment of the monomer 
density profile 



(Zm) 



/o°° Zprn{z)dz 



(5) 



Pm{z)dz ' 

in which Pm{z) is the number density of monomers as a 
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function of the distance from the grafting surface. For a 
better monitoring of the statistics of the chains conforma- 
tions, we calculate the histogram of the mean end-to-end 
distance of the chains, P{R), in which R = jj X^i^i l^-il 
and Ri is the end-to-end vector of chain i. We also 
calculate the histogram of the average distance of the 
end monomers of the chains from the grafting surface, 

P{zend), in which z^nd = jj Y^zLi and Zi is the z com- 
ponent of the end monomer of chain i. With the same 
method discussed in Ref. [1^, it has been checked that 
our results are not affected by finite-size effects (see Sec. 



B. Results 

The average thickness versus the grafting density 
for brushes of flexible, semiflexible and rodlike diblock 
polyampholytes are shown in Fig. [21 Dependence of the 
average thickness on the grafting density can be described 
well by a linear function for all values of the bending en- 
ergy, fchendj that we use in our simulations. Also, it can 
be seen that at all values of the grafting density the av- 
erage brush thickness decreases with increasing the flexi- 
bility of the chains. Linear fits to the average thickness of 
the brushes of semiflexible and rodlike chains for which 
the persistence length, Ip, exceeds their contour length, 
Lci are of approximately the same slope (see solid and 
dashed lines in Fig. The slopes of the linear depen- 
dence in cases of two brushes of flexible chains are also 
approximately the same and differ by a factor of ^ ^ from 
those of two other cases (dotted and dash-dotted lines in 
Fig. To analyze such dependencies of the brushes 
thickness on the grafting density, we look at the equilib- 
rium conformations statistics of the chains at different 
values of pa- In Figs. [3] and |4] the histograms P{R) and 
P{zend) at three different grafting densities are shown 
for brushes of flexible, semiflexible and rodlike chains. 
Because of very similar behaviors of the histograms in 
kbend — and fcbend — ^ksT cascs, the histograms of 
the brush with ktend = are not shown for clarity of the 
figures. As it can be seen in Fig. [31 in the case of the 
brush of flexible diblock polyampholytes {kbend = ^ksT), 
dependence of the histogram profile on the grafting den- 
sity is very weak. Also, in this case, the contribution of 
large values of r (r ~ Lj,) in the histogram is negligible 
which shows that the chains are mostly coiled at all graft- 
ing densities (see a sample configuration of the brush at 
Pa<J^ = 0.08 in Fig. [T^). In the case of the brush of 
semifiexible chains {kbend = '^^ksT), with increasing the 
grafting density, the values of the histogram correspond- 
ing to smaller values of r become nonzero showing that 
polyampholyte chains take buckled conformations at high 
grafting densities [2§]. As it is expected, the histogram 
P{R) of the brush of rodlike chains {kbend — 200fcsT) 
exhibit no noticeable buckling of the chains. The his- 
tograms P{zend) in Fig. [Ushow that in the brush of flex- 
ible chains at all grafting densities the end monomers are 



mostly distributed near the grafting surface showing that 
the positive blocks of the chains are mostly turned back 
towards the anchored negative blocks. The profile of this 
histogram also doesn't depend noticeably on the grafting 
density. The histograms P{zend) for brushes of semiflexi- 
ble and rodlike chains show that two maxima appear and 
their height increase with increasing the grafting density. 
By combining the information obtained from histograms 
P{R) and P{zend) for the brush of rodlike chains it can 
be understood that at high grafting densities a fraction 
of the chains which are perpendicular to the grafting sur- 
face coexist with the remaining fraction which fluctuate 
in the vicinity of the grafting surface. In the case of the 
brush of semiflexible chains also, the fraction of perpen- 
dicular chains to the brush surface coexist with the chains 
which are buckled towards the grafting surface 29] . His- 
tograms P{R) and P{zend) show that conformations of 
the brushes of semiflexible and rodlike chains change no- 
ticeably with changing the grafting density despite the 
case of the brush of flexible chains. 



III. LINEAR DEPENDENCE OF THE BRUSH 
THICKNESS ON THE GRAFTING DENSITY: 
SCALING APPROACH 

As mentioned in the previous section, our MD simula- 
tions show that the brush thickness is a linear function 
of the grafting density regardless of different values of 
grafted diblock polyampholytes bending energy. Also, as 
it is shown in Fig. [21 in cases of the brushes of flexible 
chains the slope of the fitted line to the average thickness 
versus the grafting density is considerably smaller than 
that in two other cases. We use here a simple scaling 
theory similar to that of the solution of charge-symmetric 
diblock polyampholytes |3^] to describe the linear depen- 
dence of the brush thickness on the grafting density. This 
scaling analysis is applicable to the brushes of fiexible di- 
block polyampholytes. 

Consider M flexible diblock polyampholyte chains end- 
grafted to a flat surface of area A. The degree of poly- 
merization of each chain and the fraction of charged 
monomers are N and / respectively. Electrostatic at- 
tractive interactions between oppositely charged blocks 
of the chains lead them to form a dense layer of posi- 
tive and negative monomers of average thickness h near 
the grafting surface (see Fig. [5]). Inside the layer, to 
use the blob concept we define correlation length ^ as 
a length scale that electrostatic interactions don't per- 
turb the chains statistics at smaller length scales and 
are dominant over thermal fluctuations at larger length 
scales. Let suppose that the number of monomers inside 
the correlation blob is g and we have ^ « bg'^ where h 
and V are the Kuhn length and the Flory exponent re- 
spectively. Accordingly, the layer is a melt of correlation 
blobs in which positive blobs with a high probability are 
surrounded by negative blobs. Electrostatic interaction 
between any two neighboring blobs is of the order of the 
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thermal energy, fc^T 

keT I |« keT. (6) 

Thus we obtain ^ « IsPg^ ~ H^-^b~)"^ ■ ^ result, 
the local monomer concentration inside the blob is 

Considering that correlation blobs are space filling, the 
global monomer concentration, pgiobai — = Pa^, ap- 
proximately equals to the local monomer concentration, 
Pgiobai ~ Piocal- Thus, dependence of the thickness of a 
brush of diblock polyampholytes on the grafting density 
can be obtained as 

h^p.Nh'i^^f)'^. (8) 

The main prediction of this equation in the range of its 
validity is as follows. The brush thickness, ft, is a lin- 
ear function of the grafting density, pa, irrespective of 
the values of the system parameters. Only the slope of 
this linear function, a = Nh^{^-f-p)~^^ , depends on the 
values of the system parameters. 

One should note here that despite the bulk solution 
of flexible diblock polyampholytes, in a polyampholyte 
brush the fact that the chains are end grafted to the sur- 
face introduces an additional length scale to the system, 

— 1/2 

namely d = pa .In the scaling analysis presented here, 
if the size of the correlation length, ^, exceeds the separa- 
tion between the grafting points, d, the scaling approach 
becomes inconsistent. Although linear dependence of the 
brush thickness on the grafting density is observed at all 
values of the bending energy used in our simulations, the 
condition ^ < d is valid only in the case of the brush of 
flexible chains with kbend = 0. In this case the values of | 
corresponding to lowest and highest values of the graft- 
ing density we have used in our simulations are 0.19 and 
0.45 respectively. For the set of parameters used in our 
simulations, in the case of the brush with ktend = the 
upper limit of the validity of the scaling method ~ d) 
corresponds to the grafting density pacr'^ — 0.55 which is 
quite far from our range of the grafting density. For a 
brush of flexible chains with kbend — using the values 
& ~ cr and — I gives the slope a/ ~ 35 for the brush 
thickness versus the grafting density for the used set of 
system parameters. As it is shown in Fig. [51 the value 
of q;/ obtained from our simulations of flexible chains is 
a/ ~ 20.0 which is in reasonable agreement with predic- 
tion of the scaling method. 

Our scaling method is not consistently applicable for 
the brush of flexible chains with Ip — 5cr because the 
range of the grafting density used in our simulations is 
higher than the validity range of the condition ^ < d 
in this case. This method is not also applicable for 
the brushes of semiflexible and rodlike chains because in 



these cases the Kuhn length exceeds the contour length of 
the chains. The fact that the average thickness linearly 
depends on the grafting density at all values of kbend 
used in our simulations shows that although the scaling 
method presented here can not be used to describe all the 
simulation results, this linear dependence persists over a 
wide range of the system parameters. 



IV. CONCLUSIONS AND DISCUSSION 

Brushes of flexible, semiflexible and rodlike diblock 
polyampholytes have been studied using MD simulations 
and a scaling analysis has been presented to describe the 
results of the simulation of flexible chins brush. The av- 
erage thickness as a function of the grafting density and 
histograms of equilibrium conformations of the brushes 
are obtained. Strong dependence of the system confor- 
mations on the grafting density and separation of the 
chains into two coexisting fractions at high grafting den- 
sities have been observed in cases of the brushes of semi- 
flexible and rodlike chains. In cases of the brushes of flex- 
ible chains, single-chain behavior is dominant and depen- 
dence of the brush conformations on the grafting density 
is very weak. In spite of above mentioned differences, it 
has been observed that dependence of the average brush 
thickness on the grafting density is linear for brushes of 
all different chains. This linear dependence resulted from 
our MD simulations has been described well using a sim- 
ple scaling method in the case of the brush of flexible 
chains. 

Brushes of polyelectrolytes and polyampholytes are 
dense assembly of these macromolecules in which the 
interplay between electrostatic correlations, strong ex- 
cluded volume effects and bending elasticity of the chains 
determine equilibrium properties of the system. The 
main differences between brushes of polyelectrolyte and 
polyampholyte chains originates from opposite trends 
of inter- and intra-chain electrostatic interactions and 
different rules of counterions osmotic pressure in these 
brushes. In a brush of polyelectrolyte chains most of 
counterions are contained inside the brush and their os- 
motic pressure tends to increase the brush thickness. In a 
brush of overally neutral polyampholytes however coun- 
terions are outside the brush and have no effect on the 
brush thickness. Linear dependence of polyelectrolyte 
brush thickness on the grafting density has theoretically 
been described El, 13 ■ The results of a recent sim- 
ulation of semiflexible polyampholytes ^2§\ and our sim- 
ulations and theoretical analysis here show that linear 
dependence of the average thickness on the grafting den- 
sity is also the case in the brush of polyampholyte chains. 
Amount of the flexibility of the chains which control the 
strength of excluded volume effects and inter-chain cor- 
relations, causes the equilibrium conformations of the 
brushes of three different polyampholytes to be different. 
In brushes of semiffexible and rodlike polyampholytes, 
strong dependence of the equilibrium conformations on 
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the grafting density and separation of the chains into two 
coexisting fractions at high grafting densities are resulted 
from strong electrostatic and excluded volume correla- 
tions between the chains. Similar phenomenon has been 
observed in the brush of rodlike polyelectrolytes [22j. 

To pay attention to possible finite-size effects in our 
simulations, we have calculated the average size of the 
chains lateral fluctuations as well as the histograms 
P{zend) and P{R) for a larger brush containing M = 
8 X 8 = 64 semiflexible chains at the grafting density 
Pa<7^ — 0.1 [29|. The average size of lateral fluctua- 
tions is defined as {kat) = -mYaLi{^A\)^ where Ri\\ = 
|Ri — Ri.z| and (...) denotes averaging over equilibrium 
configurations. We have found that the histograms P{R) 
and P{zend) do not change noticeably with increasing the 



system size [29j. Also, we have found that the value of 
(kat) is smaller than the lateral size, L, of the simulation 
box of a brush containing M = 25 chains at the same 
grafting density. This result shows that in simulation 
of the brush of M = 25 polyampholytes, the chains do 
not overlap with their own images. Accordingly, to avoid 
time consuming simulations of larger brushes, we con- 
centrated on the simulations of the brushes containing 
M = 25 chains. 

Acknowledgments 

We would like to acknowledge A. Naji for his useful 
suggestions. 



[1] S. Forster, M. Schmidt, Adv. Polym. Sci. 120, 51 (1995). 
[2] J. L. Barrat, J. F. Joanny, Adv. Chem. Phys. 94, 1 
(1996). 

[3] A. V. Dobrynin, R. H. Colby, M. Rubinstein, J. Polym. 

Sci. Part B: Polym. Phys. 42, 3513 (2004). 
[4] J. Riihe, M. Ballauff, M. Biesalski, P. Dziezok, F. Crohn, 

D. Johannsmann, N. Houbenov, N. Hugenberg, R. Kon- 

radi, S. Minko, M. Motornov, R. R. Netz, M. Schmidt, 

C. Seidel, M. Stamm, T. Stephan, D. Usov, H. Zhang, 

Adv. Polym. Sci. 165, 79 (2004). 
[5] Y. Hong, R. L. Legge, S. Zhang, P. Chen, Biomacro- 

molecules. 4, 1434 (2003). 
[6] S. Jun, Y. Hong, H. Imamura, B. Y. Ha, J. Bechhoefer, 

P. Chen, Biophys. J. 87, 1249 (2004). 
[7] R. Messina, Eur. Phys. J. E 22, 325 (2007). 
[8] C. M. Wijmans, E. B. Zhulina, Macromolecules 26, 7214 

(1993). 

[9] E. Lindberg, C. J. Elvingson, Chem. Phys. 114, 6343 
(2001). 

[10] J. Klos, T. Pakula, Macromolecules 37, 8145 (2004). 
[11] A. S. Almusallam, D. S. Sholl, Nanotechnology 16, S409 
(2005). 

[12] P. Pincus, Macromolecules 24, 2912 (1991). 
[13] E. B. Zhulina, O. V. Borisov, Macromolecules 29, 2618 
(1996). 

[14] E. B. Zhulina, J. K. Wolterink, O. V. Borisov, Macro- 
molecules 33, 4945 (2000). 

[15] O.V. Borisov, T.M. Birstein, E.B. Zhulina, J. Phys. II 2, 
63 (1992). 

[16] E.B. Zhulina, O.V. Borisov, J. Chem. Phys. 107, 5952 



(1997). 

[17] A. Naji, R.R. Netz, C. Seidel, Eur. Phys. J. E. 12, 223 

(2003) . 

[18] H. Ahrens, S. Forster, C.A. Helm, N.A. Kumar, A. Naji, 
R.R. Netz, C. Seidel, J. Phys. Chem. B 108, 16870 

(2004) . 

[19] A. Naji, C. Seidel, R.R. Netz, Adv. Polym. Sci. 198, 149 
(2006). 

[20] F. S. Csajka, C. Seidel, Macromolecules 33, 2728 (2000). 

[21] C. Seidel, Macromolecules 36, 2536 (2003). 

[22] H. FazU, R. Golestanian, P. L. Hansen, M. R. Kolahchi, 

Eur. Phys. Lett 73, 429 (2006). 
[23] O.V. Borisov, T.M. Birstein, E.B. Zhuhna, J. Phys. II 1, 

521 (1991). 

[24] N. P. Shusharina, P. Linse, Eur. Phys. J. E 4, 399 (2001). 
[25] N. P. Shusharina, P. Linse, Eur. Phys. J. E 6, 147 (2001). 
[26] A. Akinchina, N. P. Shusharina, P. Linse, Langmuir 20, 
10351 (2004). 

[27] A. Akinchina, P. Linse, Langmuh 23, 1465 (2007). 
[28] P. Linse, J. Chem. Phys 126, 114903 (2007). 
[29] M. Baratlo, H. FazU, Eur. Phys. J. E 29, 131 (2009). 
[30] H.J. Limbach, A. Arnold, B.A. Mann, C. Holm, Comp. 

Phys. Communications 174, 704 (2006). 
[31] G.S. Crest, K. Kremer, Phys. Rev. A 33, 3628 (1986). 
[32] R. Strebel, R. Sperb, Mol. Simul. 27, 61 (2001). 
[33] A. Arnold, C. Holm, Comput. Phys. Commun. 148, 327 

(2002). 

[34] N. P. Shusharina, E. B. Zhulina, A. V. Dobrynin, M. 
Rubinstein, Macromolecules 38, 8870 (2005). 



1.5- 



1.2- 



0.9- 



0. 



0.6- 



0.3- 



0.0- 



(a) 



k. .=200 k„T 
bend B 



• k. .=25 k„T 
bend B 



^ k. .=5 k„T 
bend B 



fiinniiir n^lifii 



^ i- *' 



10 



rio 



20 



-1 

30 



1.5 -r 



1.2- 



0.9- 



a. 



0.6- 



0.3- 



(b) 



k. .=200 k„T 
bend B 



k. .=25 k„T 
bend B 



ku j=5 k„T 
bend B 




1.5 -T 



1.2- 



0.9- 



ST 



0.6- 



0.3- 



(C) 



k. .=200 k„T 
bend B 



k. =25 k„T 
bend B 




FIG. 3: (Color online) Histogram of the average end-to-end 
distance of the chains for brushes formed by flexible, semiflex- 
ible and rodlike diblock polyampholytes at grafting densities 
a) paa^ = 0.02, b) p„cr^ = 0.06 and c) p^a^ = 0.10. 
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FIG. 4: (Color online) Histogram of the average distance 
of the end monomers of the chains from the grafting sur- 
face for the brushes of flexible, semiflexible and rodlike di- 
block polyampholytes at grafting densities a) paO"^ ~ 0.02, b) 
Paa^ = 0.06 and c) pa^^ = 0.10. 



FIG. 5: (Color online) The blob picture of a flexible diblock 
polyampholyte brush. The blobs of positive charge are more 
probably surrounded by blobs of negative charge. Positively 
and negatively charged monomers are shown by red (+) and 
blue (— ) sphere respectively. The neutral monomers are not 
shown for clarity, h and ^ are the brush thickness and corre- 
lation blob size respectively. 



